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The kagome lattice has coordination number 4, and it is mechanically isostatic when nearest neigh- 
bor (NN) sites are connected by central force springs. A lattice of N sites has 0(s/~N) zero- frequency 
floppy modes that convert to finite- frequency anomalous modes when next-nearest-neighbor (NNN) 
springs are added. We use the coherent potential approximation (CPA) to study the mode structure 
and mechanical properties of the kagome lattice in which NNN springs with spring constant k are 
added with probability V = Az/4:, where Az = z — 4 and z is the average coordination number. 
The effective medium static NNN spring constant K m scales as V for P<c and as V for V ^> n, 
yielding a frequency scale lj* ~ Az and a length scale I* ~ (Az)~ . To a very good approximation 
at at small nonzero frequency, K m ("P, uj)/hi m (V, 0) is a scaling function of lu/lj*. The Ioffe-Regel 
limit beyond which plane- wave states becomes ill-define is reached at a frequency of order ui* . 

PACS numbers: 61.43.-j, 62.20.de, 46.65.+g, 05.70.Jk 



I. INTRODUCTION 

Understanding the nature of mechanical stability, how 
it arises, and how it can be controlled is important to 
fields ranging from civil engineering to biology. James 
Clerk Maxwell was among the first to address this prob- 
lem with mathematical rigor in his 1864 paper [l| in 
which he established the conditions for the mechanical 
stability of a frame of points connected by lines or struts 
of fixed length. His approach, which can easily be gener- 
alized to treat particles interacting via central-force po- 
tentials or bonds between particles interacting via an an- 
gular potential, has since been used to study network 
dasses 043 , randomly packed spheres near jamming 
networks of semi-flexible polymers @, Q , regular 
and random [Io| pe riodic lattices, and a host of en- 
gineering problems I n this paper, we study the 
mode structure of a particular lattice - the kagome lattice 
shown in Fig. Q] - in which nearest neighbor (NN) sites 
are connected by harmonic central-force springs of spring 
constant k and next-nearest-neighbor (NNN) sites are 
randomly connected by similar springs but with a differ- 
ent spring constant re. In the absence of NNN springs, 
finite realizations of this lattice are just on the verge of 
mechanical stability in the Maxwell sense. The addition 
of NNN springs renders the system stable. When these 
springs arc added uniformly with a springconstant k that 
is allowed to approach zero continuously |9l. Il3j . the result 
is a kind of critical point at k = with a rich mode struc- 
ture and diverging length and vanishing frequency scales 
as k — > that have direct counterparts in the properties 
of randomly packed spheres near jamming (l4j . Here we 
use the coherent potential approximation or CPA [TH, [l6[ 
(otherwise know as the effective medium approximation) 
to explore in detail the static and dynamical mechani- 
cal properties and mode structure of the kagome lattice 
as a function of the probability V that NNN sites are 
connected by springs. We find that the effective medium 
NNN spring constant n m (V,uj) scales as V" 1 for small V , 
and it gives rise to length and frequency scales I* ~ 'P^ 1 



and ui* ~ V . Beyond a characteristic frequency lo* d ~ V, 
the imaginary part of a rises quickly and phonons scat- 
ter strongly. The Ioffe-Regel limit beyond which certain 
plane wave states become ill-defined is reached at a fre- 
quency slightly larger than lo* d . 



A. Mechanical Stability and Floppy Modes 

Maxwell [l[ argued that each constraint in a system 
of N points in d-dimensions reduces the number, dN, 
of zero frequency modes of free points by one. Thus, 
the number of zero frequency modes in a system with 
N c constraints is N = dN - N c . Of these, d(d + l)/2 
are those of rigid translations and rotations, leaving 
Nj = dN — N c — d(d + l)/2 internal zero frequenc 
modes, which are usually refereed to as floppy modes 
or mechanisms in the engineering literature [l2j], when 
dN - N c - d(d + l)/2 > 0. Mechanical stability re- 
quires that there be no floppy modes, or equivalently that 
dN - N c - d(d + l)/2 < 0. Though Maxwell considered 
general frames of points, his analysis can immediately be 
applied to the determination of mechanical stability of 
extended crystals and glasses - as was first done in the 
context of network glasses, which are unstable in the ab- 
sence of bending forces, by Phillips and Thorpe [U, [l?} 
- in which points are replaced by particles or atoms and 
struts are replaced by inter-particle potentials. 

If each particle in an extended system is on average 
connected by central-force potential creating bonds to 
z other particles, then, because each bond is shared by 
two particles, the system experiences zN/2 constraints 
if bond lengths are restricted to be at their equilibrium 
length and there are no redundant bonds (l8j , and it has 
Nf=dN- \zN - \d(d + 1) floppy modes. Thus, in the 
limit of large N, a central- force system is mechanically 
stable if z > z c = 2d. In an infinitely extended system, z 
is independent of N . In a finite system cut from an infi- 
nite one, particles at the boundary have fewer neighbors 
than those in the interior, and z (which is an average 
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quantity) is of order N~ x / d smaller than it would be in 
the ideal infinite system. Systems with z = z c are called 
isostatic. A finite piece of N particles cut from an in- 
finite isostatic has a surface with of order N( d ~ 1 ^ d cut 
bonds, and thus that many fewer constraints and that 
many more floppy modes. These modes generally re- 
main floppy even if they have finite amplitudes if there 
are free boundary conditions. Because of their geometry 
with straight segments extending throughout their inte- 
rior, the square, the kagome, and related lattice of N sites 
with periodic rather than free boundary conditions have 
jy(d-i)/d i nnn itesimal floppy modes - modes that have 
zero energy for infinitesimal displacements but finite en- 
ergy for large displacements [12]. The simplest example 
of an infinitesimal floppy modes is that of a mass con- 
nected to two identical co-linear springs that are rigidly 
attached to walls on either side of the mass. If the springs 
are not under tension, infinitesimal displacements of the 
particle perpendicular to the spring cost no energy but fi- 
nite displacements do. Thus the harmonic kagome lattice 
has \/^V floppy modes. 

The recently most explored random isostatic system 
is perhaps the system of randomly close-packed spheres 
at the jamming transition [l4| in the system of friction- 
less soft spheres with one-sided repulsion. This transi- 
tion occurs when the volume faction <fi of spheres just 
exceeds the critical value <p c at which they form just 
enough contacts to support a compressional load. Thus 
at (f> c , the bulk modulus B jumps discontinuously to a 
nonzero value, the number of neighbors per particle grows 
as Az = z—z c = (4>— 0c) 1 / 2 , and the shear modulus grows 
continuously from zero as G ~ ((f) — (frc) 1 ^ 2 ~ (Az) 1 . As- 
sociated with this transition are a diverging length scale 
I* ~ (Az) -1 and vanishing frequency scale to* ~ Az 
whose behavior follows from quite general "cutting" ar- 
guments @. 



B. Periodic Isostatic Lattices 

Isostaticity is not limited to random systems. There 
are a number of periodic isostatic lattices with zjvjv = z c 
NN bonds per particle. These include the hypcrcubic 
lattice in d dimensions and the kagome lattice in two- 
dimensions and its generalization to higher dimensions 
with NN bonds occupied by springs of spring constant 
k, all with z = 2d. They can be moved off isostaticity by 
introducing springs of spring constant k on NNN bonds 
either homogeneously or randomly. 

In homogeneous systems, the isostatic limit is contin- 
uously approached by allowing k — > 0. In the square- 
lattice version of this model [2], the bulk modulus is 
nonzero at k = 0, but the shear modulus vanishes as K in 
this limit. In the kagome version, both the shear and the 
bulk moduli are nonzero at K = 0. In both models, as in 
the jamming problem, there is a divergent length scale, 
I* ~ k -1 / 2 and a vanishing frequency scale to* ~ ft 1 / 2 . 

The Coherent Potential, or equivalently the effective 



medium, Approximation (CPA) [TH, [l9l - [2l| is a powerful 
tool for the calculation of properties of random systems, 
from the electronic structure of alloys [TBI . I22I ] to the elas- 
ticity of composite materials [231-25]. It has been used 
with great success in the study of rigidity percolation 
problem (l6l . [2H l27j and, more recently, in the dynam- 
ics of off-lattice systems near a rigidity threshold [28[. 
It also provides a quantitative theory of the static and 
dynamic phonon response in the nearly isostatic square 
lattice in which next-nearest-ncighbor springs (NNN) 
are randomly added to a lattice in which all nearest- 
neighbor-site connected by equivalent springs [Ioj |. In 
this lattice, the effective-medium NNN spring constant 
K m (P, w) has a zero- frequency limit, n m (V, 0), that scales 
as V 2 ~ (Az) 2 leading to lengths and frequencies that 
scale as I* ~ (Az) -1 and lu* ~ (Az) 1 as they do in the 
jamming problem. The finite-frequency spring constant 
is a scaling function of lu/lu* at small V: K m (V,Lu) — 
K m (f, 0)h(ui/u>*). As V is increased, there is a crossover 
from nonaffine response at small V to affine response with 
K m (V, 0) = Vk at larger V, much as there is in networks 
of semi-flexible poly mers as the average polymer length 
is increased [2911301 



C. Review of Results 

In this paper, we use the CPA to explore in detail the 
static and dynamic properties of the kagome lattice as 
it is moved away from isostaticity by the random ad- 
dition of NNN springs. Our approach is identical to 
that employed in our treatment of the square lattice [Io| , 
but calculations are considerably more complicated be- 
cause each unit cell has three sites rather than a single 
one. The results of our calculations follow what we ex- 
pect from our experience with the homogeneous square 
and kagome lattices and with the random square lattice. 
The bulk and shear moduli are nonzero and proportional 
to the NN spring constant k when V — >• 0. The static 
NNN effective medium spring constant K m (V,0) scales 
as V 2 leading to I* ~ (Az) -1 and lu* ~ (Az) in agree- 
ment with cutting arguments. As is the case in the square 
lattice, n rn (V,Lu)/K, m (V,0) is a basically a scaling func- 
tion of lu/lu* at small V, but with small yet important 
deviations at small lu that describe Raylcigh scattering, 
i.e., a mean- free- path that scales as lu~ 3 in 2c?. 

The outline of this paper is as follows. In Sec. |TT] we 
review the elasticity of the homogeneous nearly isostatic 
kagome lattice. In Sec. lIIII we discuss the CPA on random 
nearly isotropic kagome lattice with the NNN bonds 
randomly occupied with probability V . In Sec. [TV] we 
discuss the results of the CPA calculation, including the 
crossover of K m from V 2 to V behavior as V increase or 
k decreases, and the rapid increase of scattering at the 
characteristic frequency u* D ~ Az. 
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tion 
particle 



= Reo + u e, thus the length of bond b between 
and £' is changed into 



Rb 



(2) 



We refer to the original space in which particle £ is at 
R,£o as the reference space, and the space after applying 
the displacement field as the target space. We consider 
harmonic potentials 



Ub = ~7^(Rb — RbR) 2 , 



(3) 



where RbR is the rest length of the bond, and kb is the 
spring constant. The length Rb can be expanded for small 
displacement u as 



FIG. 1: The kagome lattice with random additional NNN 
bonds denoted by purple dashed lines. The unit cell triangle 
is marked with filled triangles. Particle 1, 2, 3 in each unit 
cell are marked in the bottom right unit cell. 




FIG. 2: The kagome lattice and its floppy modes, with the 
reference state in gray and deformed state in red. Two of its 
floppy modes are shown in this figure marked by the yellow 
ribbons. 



II. HOMOGENEOUS NEARLY ISOSTATIC 
KAGOME LATTICE AND ITS ELASTICITY 



Expansion of elastic energy in general lattice 
models 



In this section we briefly review the elastic energy in 
central-force network models, in which the elastic energy 
U can be written as a sum of the energy of each central- 
force bond 



u = J2 u b(Rb), 



(i) 



where Rb is the length of the bond and Ub is the po- 
tential energy of the bond as a function of the length. 
We consider a displacement field on the network that 
maps particle £ which is at position R^o to a new posi- 



Rb — RbO + efao ' Ufe 

+ -^5- Ub ■ i 1 ~ e &oe&o) • u fc + 0(ul ), (4) 

where R b0 = \Re>o - Rfo|, u b = u e , - u f , and e b0 = 
(Rg/Q — R£o)/|Rfo~R^o| is the unit vector pointing along 
the bond in the reference space. Thus we have 



'(R 



bO 



:U 6 



- RbR) H 

k b e b oe b o 



fb e b0 ' u 6 

RbO 



I - e &0 e 6 o , 



u 6 , (5) 



where fb = U' b {R) = kb(Rbo — RbR) is the magnitude of 
the force on the bond in the reference space. In general 
we consider the case in which the reference state is in 
mechanical equilibrium, which means that the total force 
on each particle vanishes 



h = ^ foeba = 

b{l,l') 



(G) 



where the sum ^2 b ^ e t ,\ is over all occupied bonds con- 
nected to £. However, to capture the properties of ran- 
dom networks, which often carry residual stress, the 
length of each bond is not necessarily at its rest length, 
i.e., Rbo — RbR 7^ in general. 

The change of the elastic energy from the reference 
space to the target space of the whole system is then a 
quadratic form of the displacement field 



AU 
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CbO e bO i 



which can also been written as 



ub, (7) 



(8) 



where uJJ is the component of u;, parallel to e&o and 
is the component perpendicular to e^o- 

By doing a gradient expansion on the displacement 
field, 



u 6 = R b0 e b0 fc<9 fe u(r), 



(9) 
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where u(r) is the displacement field at position r, we 
recover the elastic energy of the continuum theory 



AU = I drK ijk id k UidiUj, 



(10) 



with 



K 



ijkl 



b 



■ ^kb^bOi^bOj 



~W~ (<% ~ e b0i e b0j) 
-KfcO 



(11) 



where the summation ^ b is over bonds connecting to 
one particle, and we are using a simple lattice with one 
particle per unit cell in this illustration. The volume of 
a unit cell is denoted by vq. 



B. Elastic energy of the kagome lattice 

The kagome lattice is a lattice with three particles per 
unit cell, and we shall use the following displacement 
vector to denote the deformation of the lattice 

,xi ^,2,:c: ^,2,y) ^£,3,x^ ^£,3,y ), (12) 

where £ labels the unit cell and (1,2,3) label the parti- 
cles in the unit cell as in Fig. [2J As we have discussed, 
the elastic energy can be expanded in small displacement 
field u, and to leading order we have the quadratic form 



AU = | ■ ■ Uf, 



(13) 



which can be built from the analysis as in Eq. (JT]), but 
generalized to the case of the kagome lattice which has 
three sites per unit cell. The matrix D is called the dy- 
namical matrix of the lattice. This clastic energy can be 
written in momentum space as 

A(7=iVu n .D_ n , I .u fl ,, (14) 



— ^u q D_ q , q , u q ,, 

q,q' 



where the Fourier transform into momentum space is de- 
fined as 



u q = ^u,e- iqR - 



TV q 
q 



(15) 



where TV is the number of unit cells. The dynamical 
matrix for the homogeneous kagome lattice with all TV TV 
bonds occupied with springs of spring constant k and all 
7VTVTV bonds with springs of spring constant k is a 6 x 6 
matrix given by 

Dq,q' = TV(5 q:q ,D q (fc, «) 



D q (fc, k) — k ^ B™' q B^ Tl _ q 

mENN 



-K 



E b 

mENNN 



NNN-dNNN 
m,q m,— qi 



(16) 



where the B vectors and their derivation are given in 

App.m 



C. The homogeneous kagome lattice and its low 
energy theory 

There are six translational degrees of freedom per 
unit cell in the kagome lattice giving rise to six phonon 
branches. Of these, three are optical branches with fre- 
quencies of order \fk, two are acoustic branches with 
sound velocities of order s/k, and one is the anoma- 
lous branch. The later three branches, which determines 
the low-energy elastic theory of the kagome lattice, have 
modes in the space spanned by the following three vectors 

V! = (1/V3)(1, 0,1, 0,1,0) 

v 2 = (1/V3)(0, 1,0, 1,0,1) 

. 1 1111, 

v * = (~7!' '^'"2'^'2)' (17) 

which correspond to two translations of the whole unit 
cell in x and y directions and the rotation of the unit- 
cell triangle around its center. The low-energy theory is 
governed by the 3x3 reduced dynamical matrix obtained 
by integrating out the three high-energy optical branches, 
as shown in App. IB 11 

For small momentum, |q| < q* H = A^j2>n/k^ the re- 
duced dynamical matrix is simply diagonalizcd by longi- 
tudinal and transverse acoustic phonons (which are lin- 
ear combinations of vi and v^) with speeds of sound 
cl = v / 3fc/4 and ct = V»/4 and the rotational mode 
with a characteristic frequency loq = a/6k at q = 0. The 
bulk modulus B and the shear modulus G are related, 
respectively, to the longitudinal and transverse sound ve- 
locities through 



cl = (B + G)/q- 4 = G/q, 



(18) 



where g is the mass density, which because there are three 
atoms per unit cell, is equal to 3 in our units. Thus, 
B = 3fc/8 and G = 3fc/16. There is only weak mixing 
between the rotational modes and the acoustic phonons, 
and the system is isotropic. 

For large momentum |q| > q* H = 4y / 3K/fc, strong 
mixing between the transverse acoustic modes and the 
rotational modes occurs, and the strong anisotropy of 
the isostatic state is retrieved. The mixing is maximal 
along q x = and symmetry equivalent directions, which 
we refer to as the isostatic directions, and the resulting 
two modes are shown in Fig. |3l The anomalous branch, 
with frequency of order k, is the lower branch of the two. 
In the limit of k — 0, the lattice becomes isostatic, the 
isotropic region is squeezed to the origin, and the anoma- 
lous modes reduce to the isostatic floppy modes with zero 
frequency along q x = (YM line in Fig. [3]) and symme- 
try equivalent directions as depicted in Fig. [2] The name 
"anomalous modes" follows the nomenclature of Ref. I6f , 



5 



M K 




K r M K q r q H Q s M q 



FIG. 3: (color online) (a) Phonon dispersion along symmetry 
directions. The dotted lines are for k = and the solid lines 
are for « = 0.02. The floppy and anomalous branches are in 
red. (b) shows anomalous and shear modes along YM and in- 
dicates characteristic frequencies and wavenumbers. Frequen- 
cies Ct>o, an d ojm are defined in the text. (From Ref. 



referring to the modes developed from the floppy modes 
as the system is moved away from the isostatic point. 
For a more detailed discussion of the low energy theory 
of the elasticity of the kagome lattice, see App. IB 11 

Of particular interest is the frequency of the anomalous 
modes in the vicinity of q x = 0. The squared frequency 
of these modes can be written as 

^ 2 (q) =wa(qv) + 4& (is) 

where c x = cl = v3fc/4. The function ui\(q y ) is plotted 
in Fig. 01 It reaches a maximum value of (uig) 2 at a 
2d saddle point at q y = Qs and a local minimum value 
of (uj m ) 2 at the zone edge q y = Qm = 2n/y/3. For 
small k, Q s ~ 4(3K/2fc) 1 / 4 , (uj*) 2 ~ 3k and {uj* m ) 2 ~ 
2k. All of the characteristic frequencies ljq > ujs > 
u>m are proportional to yfn for small K. ui\(q y ) is well 
approximated between q y = Qs and q y = Qm by 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



FIG. 4: (color online) Eigenvalue of the isostatic mode along 
isostatic directions, e.g., q x = 0, for k = 5 x 10" 4 . The 
eigenvalue of the full 6x6 dynamical matrix, ui , normalized 
by «, is denoted by the red dots, and the eigenvalue of the 
3x3 reduced dynamical matrix (|B2[) is denoted by the red 
line. The blue line represent the approximation (|20|) we used 
in the asymptotic calculation in the / in CPA. 

relation function £^.^(1, 1'): 

^ e"(-q)e"(q) 
%M) = ksTY J M ^ (q ) (21) 

where fi and v label the basis defined in Eq. (|12[) of the 
6— dimensional space of u, a labels the phonon band, and 
e"(q) is the 6-dimcnsional eigenvector associated with 
mode (a, q). This correlation function is in the zero- 
frequency limit and thus independent of dissipation of 
the system. The quantities w 2 (q) are merely the eigen- 
values of the dynamical matrix with the zero-frequency 
value of the spring constant (the effective medium spring 
constant can depend on frequency as we discuss below 
in the CPA). The diverging length scale I* ~ k -1 / 2 can 
be extracted this way from the static phonon correlation 
function in experiments. 



"A ~ n pj-[QAfw| - Qs^li ~ qy(u 2 s - uj 2 m )}, (20) 

Qm — Qs 

as is evident from Fig. 2] This relation will prove useful in 
our evaluation of integrals in our CPA analysis in Scc. lIIII 

Lengths scaling as k -1 / 2 can be extracted from the 
phonon dispersion relations in various ways. One length 
is the hybridization length l* H obtained from the hy- 
bridization wavenumber q* H = ^^Jinjk = l*^ 1 sepa- 
rating the domain of predominantly transverse phonon 
behavior at low q y from the domain of predominantly 
rotation behavior at high q y . Other lengths can be 
obtained by comparing the c x q 2 term in w 2 (q) to cu 2 j 
and cj|: q* M = u M /c x = l*^ 1 = (8/v / 6)-\/^A an d 
q s = lus/c x = Ig^ 1 = Ayjn/k. An interesting property of 
uj 2 A (q y ) is that the hybridization frequency u>* H obtained 
by setting q y = q* H in the transverse phonon frequency is 
identical to lo* s '- exq^ = w# = Wg. 

One experimentally relevant quantity is the Fourier 
transform of the finite temperature static phonon cor- 



III. THE COHERENT POTENTIAL 
APPROXIMATION ON THE RANDOM NEARLY 
ISOSTATIC KAGOME LATTICE 

The CPA is a widely used method in the study of dis- 
ordered systems 0, [H, In it, a random system is 
mapped into an effective medium with no disorder that 
is described by a Green's function with a suitable self- 
energy that can capture the effect of the disorder av- 
erage of the randomness. To achieve this, one imposes 
a self-consistency constraint that the effective medium 
Green's function perturbed by the presence of single im- 
purity in the effective medium reduces to the effective 
medium Green's function when averaged over the proba- 
bility distribution of the impurity. More specifically, the 
T-matrix of this perturbation vanishes upon averaging 
over configurations that contain and do not contain the 
impurity. 

For the case of the nearly isostatic kagome lattice, the 
effective medium has all NNN bonds occupied with an 
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effective- medium spring of spring constant c m (P, u>), and 
the effective medium Green's function is identical to that 
of a homogeneous system with k = K m (V,ui). The CPA 
procedure consists of replacing one arbitrary NNN bond 
with a new bond of spring constant k s , which takes on 
the value k with probability V (bond occupied) and the 
value with probability 1 — V (bond unoccupied). This 
procedure leads to a modified dynamical matrix 



jy = d + v, 



(22) 



where 



V q , q ,(fc, k) = (k, ^"B^, (23) 



where 1 represents the arbitrary NNN bond we have 
chosen to replace into k s . This form of V follows directly 
from the calculations leading to Eq. (|16|) . It depends on 
the wavenumbers q and q' because the perturbed system 
is not translationally invariant. 

The phonon Green's function for the effective medium 

is 



G q (w) = [w 2 I-D q ] \ 



(24) 



In the perturbed system with one bond replaced, the 
Green's function becomes 



G^( W )=[ W »I- D V]^ 



(25) 



and is no longer translationally invariant. This Green's 
function can be expanded for small V 



— NSq^Gq + G q -V qiq /-G q / 

+7J E G * v i.n- G *- v «n.«i" G i' + - • • '( 26 ) 



qi 



where we have dropped the frequency u dependence 
which is the same for every G and V. This series can be 
written as 



a' — ^^q,q'G q + G q • Tq n' • G q ' , 



(27) 



where 



Tq,q< = V q , q / + — ^V q , qi - G qi -V qi , q / 

qi 

"*"jV2 Vq,qi'G qi - V qi , q2 -Gq 2 -V q2! q' 

qi,qa 



+ •••, (28) 

is the T-matrix expressed in the wavenumber basis. 

In the CPA, the effective medium spring constant n m is 
determined by requiring that the average value of G^ q , 
be equal to -ZV<5 q . q 'G q or equivalently that the disorder 
average of the T-matrix vanish: 



The evaluation of the T-matrix is simplified by the fol- 
lowing identity, 



i£v q>qi -G qi -V 
qi 



qi,q 



(30) 



(k s Km) 



2-rtNNN 



., 1 ("ST^-dNNN t i -oNNN \ -r>NNN /oi\ 

x n B i -qi ' ^qi ' B i.qi J B i -q' ( 31) 



-(k s - Kra)V q ,q'/(K m ,u), 



(32) 



where 



f{K m ,u>) = -v y 4^2 B i-q ' ' B l 



NNN 

q : 



(33) 



with Do = V3/2 the area of the unit cell in real space and 
47r 2 /uo = 8^/^/3 is the area of the first Brillouin zone in 
reciprocal space. The integral is over the first Brillouin 
zone. The Green's function G q (w) is the phonon Green's 
function in the effective medium so it depends on K m . 
Using these relations in Eq. (f2"5)l gives 



(34) 



1 + (« a - «m)/(« m ,w) 

Thus, the self-consistency equation (|2"§| requires that 

f(K m ,w)n 2 m - (1 + Kf(K m ,oj))K m +Vk = 0, (35) 

from which one can solve for the effective medium NNN 
spring constant K m for any given V and uj. The form 
of this solution at small n m depends on the behavior 
of the function f(K m ,u>) at small K m , which is in turn 
determined by the form of the anomalous mode along 
the q x = and other isostatic directions. Details of the 
calculation of f(K m ,ui) are presented in App. [Bl 

In the following discussion unless otherwise stated, we 
use reduced units with k = 1 and lattice constant a = 1, 
and thus unitless spring constants, and elastic moduli, 
and frequencies: n/k — > k, Ga 2 /k — > G, and ui/Vk — > u. 

IV. RESULTS AND DISCUSSION 
A. CPA solution at zero frequency: static response 

We first consider the case of u = 0, which characterizes 
the static response of the system. For small V, we expect 
that the effective medium spring constant n m also to be 
small and that we can, therefore, ignore the f(K m ,Lu)K m 
term in the CPA self-consistency equation (|35|) . Using 
the asymptotic small n m limit f{K m , 0) = B/^n m , where 
B = 5(1- v/273), derived in App. [Eq. (IBT31) ]. we 
obtain the equation for n m at small V 



7>T| Ks=K + (1 - 7>)T| Ks 



0. 



(29) 
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Vk = 0, 



(36) 
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which has the solution 



K m (P,0) 



-Bk + y/B 2 K 2 + AVn 



This solution has two limits: 



AV 2 if 7>« (B 2 /4)k , 
Tk if 7>> (B 2 /4)k , 



(37) 



(38) 



where A = 1/B 2 = 3(5 + 2 v / 6)/25. In the first case, 
Ky/K m 3> K m , and the solution for K rn is obtained by 
ignoring the first term in Eq. (|3l))) : in the second case, 
the opposite is true, and K m is obtained by ignoring the 
second term in this equation. In the second case, every 
NNN bond distorts in the same way under stress, and 
response is affine. In the first case n m = AT 2 <C Vk, 
indicating that the response is nonaffine with local re- 
arrangements in response to stress. Within the CPA, 
this result emerges because of the divergent elastic re- 
sponse encoded in G (and /(K m ,0)) as n m -> (See 
App. fBj) . The nonaffine regime arises when NNN springs 
are strong enough for the second term in Eq. (f3"rJ|) to 
dominate the first. As n approaches zero at fixed V , dis- 
tortions produced by the extra bond decrease and the 
nonaffine regime becomes vanishingly small. Numerical 
solutions of the CPA self-consistency equation (|3"5j) with 
the full 6x6 dynamical matrix is plotted in Fig. [SJ along 
with a comparison to the analytical solution (|37[) and the 
two asymptotic forms in Eq. (|3"5]l . 

This crossover of the effective medium spring constant 
K m between V 2 and V is different from the nonaffine- 
affine crossover in the case of random nearly isostatic 
square lattice which describes the response of the 
square lattice to shear stress, because the shear modu- 
lus (more precisely, C44) is proportional to the effective 
medium spring constant K m for the NNN bonds in the 
square lattice. In the case of kagome lattice in this Pa- 
per, the shear modulus is finite and proportional to k, 
whereas K m mainly determines the rigidity of the lattice 
with respect to the floppy mode, which are essentially 
rotations of the triangles as shown in Fig. [2J 

Length and frequency scales can be extracted in the 
static limit much as they were extracted in the homoge- 
neous case discussed in Sec. [TTJ The finite temperature 
static phonon correlation function is the inverse of the 
dynamical matrix evaluated at w = 0, whose eigenval- 
ues and eigenvectors are identical to those of the homo- 
geneous case with k replaced by n m = K m (V,0). The 
eigenvalues allows us to identify frequencies by taking 
the square roots of the appropriate eigenvalues of D: 



J o 



= \/<OK m >Wg=: \f^K m > W 



M 



(39) 



Unlike the situation homogeneous lattices, these frequen- 
cies are not equal to any physical dynamical-mode fre- 
quency of the system. They do, however, provide infor- 
mation about the static properties of the phonon cor- 
relation function 5f that could in principle be measured 
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FIG. 5: (color online) CPA solution at zero frequency. Data 
points show the numerical solution n m as a function of V at 
uj — of the CPA self-consistency equation (|35|) with the 
full 6x6 dynamical matrix. NNN bond spring constant 
k — 1(P 2 , 10°, 10 2 are shown in red dots, purple squares, and 



blue diamonds respectively. The analytical solution (|37[) of 
these three cases at small V are shown in red, purple, and 
blue lines. Also shown are the nonaffine (n m = AV 2 ) and 
affine (n m = Vk at k = 10™ 2 ) limits in green dashed line 
and red dash-dotted line. At large V the numerical solution, 
especially the one for k = 10 2 deviate significantly from the 
nonaffine limit form because Eq. (|36|) is an approximation at 
small V by ignoring the highest order term in Eq. (|35[) . 



at finite temperature via scattering or particle tracking 
experiment. They also lead to diverging lengths just as 
they do in the homogeneous case: 



I* = I 



H 



^ = ^' = 1757 = ^ (40) 



B. 



CPA solution at finite frequency: dynamic 
response and damping 



For finite frequency uj, the effective medium spring 
constant is complex, K m (V,oj) = K' m {V,uj) — ik'^V ,oS), 
where the imaginary part K m {V,uj), which describes 
damping of phonons in this random network, is odd in 
uj and positive for uj > 0. From the analysis for the 
static limit uj = 0, we see that the interesting case is 
the nonaffine regime with V -C (B 2 /4)/c, in which the 
self-consistency equation (|36p simplifies to 



f{K m ,Uj)K m = V. 



(41) 



— 1/2 

In the static limit, f(K m ,0) ~ K m is singular in the 
K m — > limit. As we show in App. IB1 at finite frequency, 
f{n m ,uj) - [^(3k to -uj 2 )/n m - y/ '{2K m - w 2 )//c m ]//c m , 
which leads to 



K m {V,uj) = ^-( 5 + 2 ^6 



25uj 2 
18V 2 



(42) 



as depicted in Fig. [51 Taking ui = 0, this solution re- 
duces to the zero-frequency solution [Eq. (|38[) ] in the 
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nonaffinc limit. Equation (|42j) develops an imaginary 
part when \ui\ > 2\f%V /5, which, as we discussed above 
must be negative for lo > 0. It is straightforward to see 
that this solution satisfies the scaling form n rn (V,Lo) = 
K m (V, 0)h(oj/uj*), as does the CPA effective NNN spring 
constant in the square lattice Ref. This solution 

shows a rapid increase of damping beyond a characteris- 
tic frequency 



(43) 



marking another characteristic frequency that scales also 
as V. 

Numerical solution of the CPA self-consistency equa- 
tion ([33)) using the full 6x6 dynamical matrix is also 
shown in Fig. [6] We see that the asymptotic form (|42]) 
captures the solution fairly well. 

This special behavior of the imaginary part of the effec- 
tive medium spring constant K m is related to the phonon 
spectrum of the kagome lattice. As we have discussed 
in Sec HH at low frequencies, there is only very weak 
mixing between the rotational branch which is strongly 
affected by the NNN bonds and the acoustic phonon 
branches which arc only very weakly scattered by the 
NNN bonds, and thus the damping to the acoustic 
phonons is very weak. On the other hand, for higher fre- 
quencies the transverse phonon strongly mixes with the 
rotational modes, making the damping rapidly increase 
with scatterings from the NNN bonds beyond lo* d , al- 
though the exact value of lo* d may be an artifact of the 
CPA method. The weak scattering below lo* d is not cap- 
tured by the asymptotic form (1421) for small K m because 
the function n m (V,io) in Eq. (|42p was obtained using the 
dominant small n m limit of the integral f(K m ,Lo). There 
are, however, contributions to this integral that do not di- 
verge and that contribute a subdominant imaginary part 
to K m , even when lo < lo* d , that is of order V s lo 2 at small 
lo corresponding to Raylcigh scattering. More discussion 
is included in App. [B] 

In the homogeneous case, the eigenvalues of the dy- 
namical matrix lead naturally to the identification of 
characteristic frequencies u>* s , and ui* M that vanish as y/~K 
in the limit of k — > 0. In the random case, we have to 
deal with both the frequency-dependence K m (V,Lo) and 
the fact that it is a complex number, and we must ask 
whether these frequencies have any real meaning. As 
discussed earlier in Sec. Ill CI we can extract frequencies 
from the static dynamical matrix in exactly the same way 
that we did for the homogeneous case, and they satisfy 



LO e 



3.85w|5 > ui'g = 2.72a; * D > lo 



2.22lu* d > lj* d . 



(44) 

Thus all of these frequencies are greater than the fre- 
quency Lj£>. As a result, the signatures in the phonon 
dispersion relation including the hybridization and the 
saddle point are washed out by the strong scattering, as 
is shown in Fig. [5] 




FIG. 6: (color online) CPA solution at finite frequency for 
V = 0.01 and V = 0.05. The numerical solution to Eq. (|35)l 
with the full 6x6 dynamical matrix is shown as the data 
points. Blue circles and red squares represent real and (neg- 
ative of) imaginary parts of K m at V = 0.01, and green dia- 
monds and brown triangles represent the real and (negative 
of) imaginary parts of K m at V = 0.05. The asymptotic 
form (|42p is shown as the blue (real) and red (negative of 
imaginary) lines. In this plot frequency is rescaled by V, and 
the effective medium spring constant n m is rescaled by its 
value at zero frequency which is real. 



C. Phonon density of states 

The phonon density of states (DOS) can be calculated 
from the retarded Green's function through 



p(io) = TrImG(q,w) 

7T 



(45) 



where the trace is over both momentum q and the phonon 
modes. Using this we get the phonon DOS of the effective 
medium solved from the CPA, as plotted in Fig. [7] As a 
comparison, we also show the phonon DOS of a periodic 
kagome lattice with the NNN spring constant equal to 
KmX'PjO), which is real valued. 

For small frequencies, at which the imaginary part of 
the CPA solution K m ("P,u;) is very small, the two DOS 
are very close, and can be fitted nicely by the Debye-like 
total DOS of the transverse and the longitudinal phonons 



Ps{u) = 



(4tt/V3)c| (4tt/V3)c 



(46) 



where c 2 L 



3fc/16 and c\ = fc/16 are respectively the 
longitudinal and transverse speed of sound (we have 
taken k = 1 as stated earlier). 

At the critical frequency lo* d , the imaginary part of 
K m (V,io) rapidly increase, inducing a rapid increase of 
the phonon DOS. On the other ha nd, the per iodic lattice 
exhibit a jump in DOS at lo* a[ = y / 2K m (7- > , 0) correspond- 
ing to the minimum of the phonon dispersion relations 
at the edge of the 1BZ. At w* = yJZK m (V,Q) the DOS 
of the periodic lattice has a logarithmic singularity, cor- 
responding to the saddle point of the phonon dispersion 
at Qs — 4(3K I „/2fc) 1 / 4 on the isostatic directions 0. 
For the CPA effective medium, this singularity is totally 
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FIG. 7: (color online) (a) The phonon DOS at V = 0.01 (blue 
circles) and V = 0.05 (green squares) of the CPA effective 
medium and the pure kagome lattice with the NNN spring 
constant equal to the zero frequency effective medium value 
K m (V, ui) for V = 0.01 (purple line). The frequency is rescaled 
by V . (b) The phonon DOS at small frequency for CPA ef- 
fective medium (color scheme the same as in (a)). The Debye 
DOS defined in Eq. (|46[) is also shown as the red dashed line. 

washed out due to the strong damping beyond ui* D . which 
is similar to the case of the square lattice fioj . 

D. Phonon scattering and the Ioffe-Regel limit 

From the CPA solution at finite frequency, we identi- 
fied a frequency scale ui* D beyond which phonon scatter- 
ing rapidly increase. In this subsection we examine the 
scattering of phonons in more detail. 

The scattering of the transverse phonons is character- 
ized by the imaginary part of the phonon response func- 
tion projected to the transverse direction Im%rr(q, a;). 
The phonon response function is defined as 

= Iffy w 

where t and t' label time, [i and v label the basis de- 
fined in Eq. (fT2|) of the 6— dimensional space of u. This 
response function is related to the phonon Green's func- 
tion through x = —G. The imaginary part of the trans- 
verse component of this response function ImxTTiQ, 




FIG. 8: (color online) Scattering of phonons characterized 
using the imaginary part of the transverse component of the 
phonon response function ImxTr(q, i^) as a function of ui for 
various values of q y (we took q x = to follow the isostatic 
TM direction). The green line in the bottom plane marks 
LUp, the blue and red lines marks the solution u/ and uj" of 
the equation uj 2 — uja(<1, K m (V,oj)) 2 = 0, which solves for the 
pole of the Green's function for the anomalous branch. The 
derivation of uja is shown in App. [B"T1 



characterizes the scattering of the transverse phonon by 
disorder. IrnxTT^, uj) is calculated for small momentum 
and frequency using the asymptotic CPA solution (|42p . 
and shown in Fig. [5J Also shown in the figure is the 
frequency at which the phonon Green's function of the 
anomalous branch has a complex pole, which is solved 
from the equation ui 2 — uja (q,K m (P,uj)) 2 = 0, which 
characterizes the dynamic dispersion relation. We use 
the form of uja as defined in App. IB II for this calcula- 
tion. It is clearly shown in the figure that below uj* Di the 
response function has Dirac-delta peaks at the frequen- 
cies determined by the transverse phonon dispersion re- 
lation uj = CTQy Above ui* D the imaginary part uj" takes 
off, and the phonon peaks progressively broaden, show- 
ing that the transverse phonon is no longer a good eigen- 
state of the system. Furthermore, the hybridization of 
the transverse phonon and the rotational modes into the 
anomalous mode and the van Hove singularities in the 
density of states are washed out by the strong scattering. 
As a result, ujq, ujg, ui* M no longer play meaningful role 
in the dynamic response function. 

The strength of the scattering can be characterized by 
the Ioffc-Rcgcl (IR) limit, which states that the plane 
wave states are no longer well defined if the mean free 
path i m f p is comparable to or less than the phonon wave- 
length A. An equivalent condition is that the relaxation 
time becomes comparable to the period of the wave, i.e., 
bj" ~ uj' . The solution for the positions of the complex 
poles of the Green's function of the anomalous branch 
uj 2 — ^(q, K m (V,ui)) 2 = shows that the imaginary 
part to" becomes comparable to the real part ui' not far 
beyond ui* D and that Wj R ~ ui* D ~ Az. The associated 
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IR length scale can be derived from uj* d and ct to be of 
order ~ y/k/n m ~ Az . 

This IR length scale differs from that in j ammed solids, 
Id ~ Az -1 / 2 as derived in Rcfs. [3l|, [32| (called l s in 
Ref. [28| ) . This discrepancy can be attributed to the dif- 
ferent scaling of the shear modulus G. In the kagome 
lattice, the shear modulus G is proportional to k and 
thus scales as Az°, whereas in jammed solids G ~ Az. 
Thus, the transverse speeds of sound scales as (Az)° 
and (Az) 1 / 2 in these two cases respectively. In both 
cases, the frequency beyond which plane wave states are 
strongly scattered is a;* ~ Az. Therefore, the scatter- 
ing length scale, ct/uj* are respectively I* ~ (Az) -1 and 
Id ~ (Az) -1 / 2 in the kagome lattice and jammed solids. 



E. Comparison between different random nearly 
isostatic systems 

Up to now, three examples of random nearly isostatic 
systems have been studied, including the random nearly 
isostatic square lattice, random nearly isostatic kagome 
lattice discussed in this Paper, and jammed solids near 
point J. In all cases, the characteristic frequency for the 
onset of the anomalous mode plateau uj* ~ Az and the 
isostatic length scale I* ~ (Az) -1 . On the other hand, 
the scaling of clastic moduli in the three cases are dif- 
ferent because of different network architecture: in the 
square lattice G ~ n m ~ (Az) 2 and B ~ k ~ (Az)°, 
in kagome lattice G, B ~ k ~ (Az)°, and in jammed 
solids G ~ Az and B ~ ~ (Az)° (a jump from 
zero to finite value at point J). As a result, the scat- 
tering length scales arc also different. The square lattice 
is anisotropic, and we studied the scattering of the u x 
vibrations along q y direction and found that the scat- 
tering length corresponded to the point M in the first 
Brillouin zone, q y = n, or l x ,m. = «• In kagome lattice 
Iir ~ ct/uj* ~ (Az) -1 . In jammed solids, the IR length 
scale Id ~ ct/uj* ~ (Az) -1 / 2 . 

In Ref. [23], Wyart studied transport properties of 
amorphous solids modeled by an isostropic random net- 
work near its percolative rigidity threshold. In this sys- 
tem, both B and G vanish as Az, and the crossover fre- 
quency between plane- wave and strongly scattered states 
is a;* ~ Az. Both the longitudinal and transverse sound 
velocities scale as (Az) 1 / 2 , and the IR length Jtr ~ 
cl.t/uj* scales as (Az) -1 / 2 for both modes. The CPA 
self-consistency equation for low frequency in Ref. [28j . 
Eq. (7), can be rewritten in the form 

k 2 M -(Az)k m +u 2 = 0, (48) 

by ignoring the w 3 term which correspond to Rayleigh 
scattering in the second subequation. The solution to 
this equation 




has a form very similar to the that of the effective medium 
NNN spring constant n m in the kagome lattice as shown 
in Eq. (j42|) . Although fejwr (c*j = 0) ~ Az in the amor- 
phous solid and kai(uj = 0) ~ (Az) 2 in the kagome lat- 
tice scale differently with Az, the frequency dependence 
of kM(w)/kM(oj = 0) and km(uj)/km(uj = 0) are almost 
identical: they are both of the form a + bWl — (oj/uj* d ) 2 
where a and b are constants and Ujj ~ Az. The w 3 term 
ignored in the above analysis leads to Rayleigh scattering 
and corresponds to the subdominant terms in / discussed 
in Sec. IIVBI The difference between w 3 in Ref. [28[ and 
w 2 for the kagome lattice is the spatial dimension. 

To summarize, we examined the random nearly iso- 
static kagome lattice via the CPA, we obtained effective- 
medium NNN spring constant K m that scales with the 
occupancy probability V ~ Az of the NNN bonds as V 2 
at small V. Below the characteristic frequency lo* d ~ V, 
there is only weak damping of acoustic phonons arising 
from Rayleigh scattering, whereas above uj* d scattering 
increases rapidly and the system shows proximity to the 
IR limit. We compare the kagome lattice to other nearly 
isostatic systems including the square lattice, jammed 
solids near point J, and a model random isotropic net- 
work [l^. The characteristic frequency scale u>* ~ Az, 
marking both the onset of the plateau of the anomalous 
modes and the strong scattering of plain wave states, is 
found to be a universal property of all of these systems. 
The elastic modulus G, B and thus the transport length 
scale depends on the network architecture and are not 
universal. 
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Appendix A: The dynamical matrix of the kagome 
lattice 

To construct the dynamical matrix of the kagome lat- 
tice, we use the form of the elastic energy given in Eq. (J7J. 
Because we consider the reference state of all bonds at 
their rest length, we have /& = 0, thus there is only pro- 
jection of u onto the direction along the bond. We first 
consider the case of simple lattice with one particle in 
each unit cell and rewrite Eq. (J7J) as 

b 

b 

x(<Wi - &',< 2 ) e <i<a ' U£f, (Al) 

where £1,(2 labels the two particles connected by the 
bond b. Thus the dynamical matrix D, as defined in 
Eq. (|13|). is given by 

b 

x(^',fi - 5e>,e 2 )e ei £ 2 . (A2) 
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Dqq' 



It is convenient to express the dynamical matrix in 
momentum space via the Fourier transform defined in 
Eq. CS) 



i 

= E e_iq ' r ' w ' r{ 'EE^ e ^^ 

x {fie,ti — 8£,£ 2 )(fii',ei ~ $t' ,e 2 ) e £ie 2 

b 

x(l- e ^' b )e b e b 



(A3) 



where the ' above the summation of £2 denote a summa- 
tion over particles connected to £\ , and b = r> — rep- 
resent the bonds connected to an arbitrary particle (note 
the difference from b in the previous equation, which rep- 
resent all bonds in the system). One can define the dy- 
namical matrix for translational invariant system as 



D q,q' 
D„ 



E^b(l-e- lq ' b )(l-e^' b )e b e b 

b 

^ ^ & m B m ,qB m ,_q, (A4) 



where the summation m is over bonds connected to an 
arbitrary particle, and the vector 



B m , q = (l-e-^ b '")e b 



(A5) 



is a convenient way to express the dynamical matrix. 

For the kagome lattice, which has three particles per 
unit cell, one need to modify the above construction of 
the dynamical matrix, and in the basis of 

u £ = (u^.i j2 ;, U£ t i : y, U£.2,xj u l.2,y, u t,3,x: u £,3,y ), (A6) 



bond is counted once) 

1 1 a/3 



NN 

q 



B 2,q 



B 3,q 



■dNN 

B 4,q 



t>N N 

B 5,q 



BNN 
6,q 



2' 2 ' 2' 2 
0,0,1,0,-1,0^ 



,0,0 



- -— 0-- —) 

2' 2 ' ' ' 2' 2 J 

2' 2 '2 

^ e -*U«-+^9») }0> 



0,0,-e-^, 0,1,0 




^ « ,-I e - i (-i«-+^9«), 



2' 2 
V3 _ 
— 6 



2 



(A8) 



and the B vectors for NNN bonds for each unit cell 



Bi 



NNN 

q 



-ryNNN 

B 2,q 



r>NNN 

B 3,q 



ryNNN 

B 4,q 



BNNN 

5,q 



B 6,q 



V 2 2 ' ' ' 2 2/ 
(o,0,0, e -*(^ + ^),0,-l) 

i _^ e -*(*-+^) 

V 2 ' 2' ' ' 2 

_L-i(l&+f%) 
2 



2 ' 2' 2 



— <iy J 



,o,o 



(o,0,0,l,0,-e^(-^+^)) 

_^ e <fa L*fa 2^ I 0> o'). 

2 2 '2 2 ;v; 



Appendix B: Calculation of the asymptotic form of 
the f(n m ,u)) function at small K m 



with particles 1,2,3 labeled as in Fig. [2J the dynamical 
matrix can be expressed as 

D q , q < = iV£ q , q 'D q (/c, k) 

D q (fe, «) = k E] B m,q B m,-q 

m£NN 

+n E B ™B™, (A7) 

meNNN 

with the B vectors for NN bonds for each unit cell (each 



1. The reduced dynamical matrix 

To calculate the asymptotic form of /(K m ,w) we first 
simplify the problem by reduce the dynamical matrix into 
the space of its three low energy modes by integrating out 
its three high energy modes [9[ . The resulting low-energy 
dynamical matrix is conveniently represented in the basis 
of longitudinal and transverse phonons and the rotational 
mode (mode ^3) 

(1/1,1/3,^3) = ( j^j . |^j ( B1 ) 
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in which the dynamical matrix takes the form 



/ 



321 
16 





2_ 

, 16 



|g cos 30 
-fjt sin 36» 



\ feces 30 -fesinSfl fg + fi^ 



in leading order of small « and quadratic order in q (the 
cross term of order K m q 2 is considered higher order and 
has been dropped). 

Eigenmodcs of the dynamical matrix arc identified by 
diagonalizing D^. Strong mixing between the trans- 
verse mode and the rotational mode occurs along q x = 
(i.e., 9 = 0) and symmetry equivalent isostatic directions. 
The resulting two eigenvalues (by diagonalizing the lower 
right 2x2 block) are 



obtained from the quadratic order of the renormalized 
3x3 matrix t)( R \ The lower eigenvalue Cj\ correspond 
to the anomalous mode, which is close to the trans- 
(B2) verse mode v' 2 (which is simply v\ for q x = direc- 
tion) for q y q* H = 4y/3n m /k . For q y > q* H this 
anomalous mode corresponds to the linear combination 
of (v' 2 — v' z )/2, which is actually the floppy mode of the 
kagome lattice in the K m — > limit, in which uia — > 0. 



h +3Kr 

16 + ^ 



(|)V(3^ 



+(3 Km ) 2 , (B3) 



J 



2. Leading order divergence of f(K m ,cj) 



The function f(K m ,u>), as given in Eq. (|33|) . can be analyzed using the simplified dynamical matrix (|B2j) . which is 
the leading order form in small K m and q. Thus we can obtain an asymptotic analytical calculation of the integral / 
by projecting from the 6-dimensional basis in Eq. (fl~2]) onto the 3-dimensional basis in Eq. (|B1[) built from the three 
low energy modes of the system 



f(n m ,uj) = - 



d 2 q 

bz 8^773 

1BZ 87r 2 /\/3 



B™-G q (o;).B™ 

B NNN . qTq . Gq(w) . qTq . B 7VA 



NNN 



where 
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<?x 1v qx 1y gx <h 

\fZq \[Zq \fZq \/Zq \fZq v3<? 

qy 9x 1y gx qy gx 

~73 2V3 ~2 2 



(B4) 



(B5) 



is the orthogonal transformation from the basis of U£ = (v,£ i x , u^i ug^.x-, u i,2,y, u e,3,x: u i,3,y) to the basis (i^, z^ 2 , ^3) 
in Eq. (|B1[) with the longitudinal, transverse, and the rotational mode. In these new basis, the dynamical matrix is 
modified by integrating out the high energy modes and keeping to leading order in small n m and which lead to 
the simple form of Eq. (|B2[) and thus the Green's function can be analyzed correspondingly. Note that we use 
the Green's function G q (w) calculated from the renormalized dynamical matrix (IB2I) . so it is different from the bare 



(B6) 



value ■ G q (u) • . The transformed the B^ ^ vector in the basis of (v[, v^v'^) takes the form 



' (e-*- - 1)(3& + y/3q y ) (e"^ - l)(>/3<fe - 3g„) 



6|q| 



6|q| 



'-2 (1 + C 



The leading order term of this integral in small n m is from 1/3, the anomalous mode, which has a small frequency of 
order ^/K m over the whole range of momentum from qu to the edge of the Brillouin zone along the isostatic directions, 
and thus correspond to diverging contributions to the / integral in small K m . 

For an approximation of the / integral at small K m , we use the dynamical matrix of the form (|B2[) . which kept to 
leading order in n m and quadratic order in q. At small momentum, the dynamical matrix (|B2[) is diagonalized by the 
basis (v[, v' 2) v' 3 ), and the Green's function G q (w) takes the form of a diagonal matrix 



1 



1 



13 



which is isotropic and valid for small momentum |q| < q* H . Thus, the small momentum region contribute to / the 
following terms 

dg x dq y j- (1 - cosq x ){3q x + \/3q y ) 2 (1 - cos q x ){ \/3q x - 3q y ) 2 1 + cos q x 

M<qll 8nyV3l 18{q 2 + q 2 ){u) 2_^+^l ) l8{q 2 + q 2 )( u2-«l±&) 2 (^-6 K -^) 



At large momentum, the dynamical matrix can be diagonalized to leading order in n m in the basis (v[, " 2 ^ 3 , " 2 ^ 3 ) i 



in which Bi q takes the form 

/ (e-<«- - l)(3q x + V3q v ) {er** - l)(VSq x - 3q v ) 1 



6|q| 672|q| 2V2 

(e~^ - l)(V3q x - 3q y ) , 1 



(1 



6V2|q| 2^2 
and the Green's function G q (w) takes the form of a diagonal matrix 



(1 + e-**-)), (B9) 



1 1 

,2 



G q (cj) = Diag 2J _ 2 , ix-rj ; 3"v)j ( B1 °) 



which is for the direction of ^ = 0, and we have used the approximated form [20] of that represent the dispersion 
relation of the anomalous mode at large frequency, as depicted in Fig. 21 For this calculation we use the small n m 
values (wj) 2 = 3« m and {oj* m ) 2 = 3n m . 

For the other two directions one should change the third term above from q 2 into the perpendicular direction of 
the two isostatic directions accordingly. Thus we need to divide the first Brillouin zone into 3 parts: \9 — tt/2\ < 7r/6, 
\6 — 7r/6| < 7r/6, and \6 — 57r/6| < 7r/6, and integrate each of them out separately and then calculate the sum. Here 
we just do the \6 — 7r/2| < 7r/6 part as an example, which uses the form of the Green's function in Eq. (IB10|) . This 
part of the integral is 

, , v 2 f^TT f (l-cosq x ){3q x + V3q y ) 2 

/>, f (K m ,w) = - — — 7= I dq y dq x \ 3 ( g 2 +0 2^ 

tofi/yfiJa J-l$ ^i 8((Z 2 + (? 2 )(w 2_H£^)) 

(2 + cos q x )ql + \/3(l - cos q x )q x q y + 3q 2 y 

(2 + cos q x )q 2 x + V3(l - cos q x )q x q y + 3gg ^ 



&(q 2 x +q 2 y ){^ 



Qm-Qs 16 



and the integral for the other two directions can be calculated similarly. 

The leading order contribution to f{K m ,u>) in small n m is from the third term in Eq. (|B11[) . which represent the 
isostatic mode. We first consider the uj = case, for which the leading order term of f(K m ,u>) is 

2 [ik {2 + cos q x )ql + V3(1- cos q x )q x q y + 3q 2 

/f,i.o.(«m,0) = - - ; 7= / dq y dq x 



8lT 2 /\/3j n * V J-\2y± X R( n 2 I „2W Q» IJ |-Qs^f-h(-'|-" , M) 1 3g'\ 
' v 1h v \1x^1y)\ Qm-Qs ^ 16 ) 

2 /"^. {2 + cos q x )ql + V3(1- cos q x )q x q y + 3ql tt 4 

87r 2 /V3 A- 9y 6( g 2+ 9 2) / QAf ^-q s ^ f - gB (a,|-^ f ) 7! 

V Qm-Qs 



2(1 - 7273) 



(B12) 



where we took the limit of K m — > and make use of the identity lim a _>.o a 'i+ x s — (^/a,)5(x). Adding up the contribution 
from 9 = it/ 6 and 6 = 5n/6 part we have 

5(1 - ^2/3) 

/i.o.(«m,0) ^ 7^= ■ (B13) 



14 



Other terms in Eq. (|B8|) and (jBllI) contribute higher order terms in small K m , and is discussed in Sec. IB 31 
In the case of w > 0, the leading order term can be calculated in a similar way 



/s,i.o.(«rn>°) = - - „ , 7= / dq y dq x 



2 


J 1h J 


8tt 2 /V3 
2 / 


8tt 2 /V3 y, 

2 (j 


»--- 


V^m \ V 



5s . (2 + cosq x )q 2 x + y/3(l - cosq x )q x q y + 3q 2 



o(q 2 + q 2 )^ 2 _ M) 

27ri(16/3) 



2(4/V3)yC 2_ ^3 



Qa/-Qs 



(B14) 



The q x integral can either be evaluated using the S function trick by assuming an infinitesimal imaginary part of u) 
(u — > lo + iS), or by extending the integral limit of q x to (—00,00) (because the integrand decays fast when q x is large) 

and using contour integral. We also assumed that |^/ Q M "s ^g"j"_q^" s " A — — uj 2 \ <C 1 to make the simplification 
that cosg x ~ 1. Adding up the contribution from 9 = 7r/6 and 9 — part we have 




/z. .(K m , W )~ -J==\ 

Other terms in Eq. (|B8|) and (jBllI) contribute higher order terms in small K m , and is discussed in Sec. 



3. Correction at small frequencies 



To get the correction to the asymptotic solution of K m (V,uj) as in (|B15|) . in particular the small imaginary part 
rather than zero at small frequency, we calculate the imaginary part of / at small frequencies and solve for the 
correction to K m (V,u) perturbatively in the CPA equation. 

Because we consider small frequencies w 2 < K m , the contribution is from the two acoustic modes, which are isotropic, 
and thus can be calculated as 



Im/i 



1h 



8ir 2 /V3 Jo Jo 



dq / d9qlm 



q 2 cos 2 6»(3cos6» + v^sinfl) 2 



18(w 2 - f^q 2 +l5) 



20 2 
27"' 



and 



Im/j 



'Ih 



&ir 2 /V3 Jo Jo 



2tv 



dq d9 q Im 



q 2 cos 2 6»(V3cos6»-3sin6l) 2 



lS( w 2 -±q 2 + i5) 



Thus we have the correction to Eq. (|F315|) that is valid for small w as 

5 



/(tm,w) 



,3 -i--\r-^ 



128 2 
+ 1 — or. 
27 



(B16) 



(B17) 



(B18) 



We then solve the leading order CPA equation in small K m nonaffme regime perturbatively using this corrected 
form of / at small uj, and get 



(o) 



256 



135(1 - v/2/3) 



(0h3/ 2 ;,„2 



(^ } ) 3/ 



(B19) 



where is the zeroth order solution (|42[) . This correction is very small and can not be observed in our numerical 
solutions within precision. 
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